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We present a few-parameter ansatz for pulses to implement a broad set of simultaneous single¬ 
qubit rotations in frequency-crowded multi-level systems. Specifically, we consider a system of two 
qutrits whose working and leakage transitions suffer from spectral crowding (detuned by 5). In order 
to achieve precise controllability, we make use of two driving fields (each having two quadratures) 
at two different tones to simultaneously apply arbitrary combinations of rotations about axes in 
the X-Y plane to both qubits. Expanding the waveforms in terms of Hanning windows, we show 
how analytic pulses containing smooth and composite-pulse features can easily achieve gate errors 
< 10“^ and considerably outperform known adiabatic techniques. Moreover, we find a generalization 
of the WahWah method by Schutjens et al. [Phys. Rev. A 88, 052330 (2013)] that allows precise 
separate single-qubit rotations for all gate times beyond a quantum speed limit. We find in all 
cases a quantum speed limit slightly below 27r/5 for the gate time and show that our pulses are 
robust against variations in system parameters and filtering due to transfer functions, making them 
suitable for experimental implementations. 


I. INTRODUCTION 

Quantum information processing and offshoots thereof, 
such as quantum computing and control are rapidly 
evolving fields in physics. Often, these involve physi¬ 
cal systems for which a computational subspace can be 
singled out of a larger Hilbert space that is feasible for 
logical operations. Candidates for systems with such a 
structure are for instance several types of superconduct¬ 
ing qubits [1], optical lattices [2], quantum dots [3, 4], 
Rydberg atoms [5, 6], neutral atoms [7], diamond NV 
centers [8], trapped ions [9, 10] and nuclear magnetic res¬ 
onance [11]. A qubit subspace in these systems is typi¬ 
cally formed by two levels of a multi-level system, such as 
an anharmonic ladder. This directly connects to leakage 
out of the computational subspace - especially on short 
timescales - since the transition frequencies of successive 
levels may sometimes only differ by a few percent [12]. 

We focus our attention on superconducting qubits em¬ 
bedded in circuit QED (cQED) architecture, where they 
are coupled to microwave cavities [13-15], both because 
it is a leading candidate for quantum computation, and 
because of the critical role that leakage and frequency¬ 
crowding play in such systems. Like the other candi¬ 
date architectures, decoherence will set a timescale below 
which gates must be performed. Decoherence may for ex¬ 
ample arise from the electromagnetic environment [16] as 
well as from intrinsic material properties [17]. Optimized 
designs of superconducting qubits, such as the transmon 
[16] and its 3D version [18, 19] avoid decoherence and al¬ 
low for coherence times approaching 100/iS, thus an im¬ 
provement by several orders of magnitude compared to 
previous designs having coherence times of a few nanosec¬ 
onds [15]. 

Experimental progress in implementing fast gates has 
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been significant, achieving fidelities exceeding 99% for 
both single-qubit and two-qubit operations [20, 21], and 
fidelities ~ 85% for two-qubit algorithms like Grover’s 
search algorithm [22]. Leakage was removed in these ex¬ 
periments using the derivative removal by adiabatic gate 
(DRAG) [23-25] technique, which works by considering 
the adiabatic effect of driving the leakage transition on 
the qubit’s transition. 

Apart from decoherence and leakage, addressability of 
single qubits in a multi-qubit layout [26, 27] is another 
key challenge for large scale quantum computing. Often 
enough it is impossible to have a single control field for 
each qubit, for instance if more than one qubit are in 
the same cavity, hence being spatially too close to be ad¬ 
dressed individually [25]. Then, the applied field collec¬ 
tively drives all qubits at once, necessitating addressing 
qubit control by an internal parameter such as carrier 
frequency. A phenomenon that arises in such cases is 
frequency-crowding [28-30], e.g. describing a situation 
where logical transitions are well separated but for in¬ 
stance the working transition of one qubit is too close to 
a leakage transition of another qubit to be driven individ¬ 
ually. In Ref. [28] a situation with two 3D transmons in 
one cavity was studied with respect to driving a rotation 
on one qubit, leaving alone the other one. They showed 
that it is possible to achieve a high fidelity by supple¬ 
menting a Gaussian pulse with sideband modulation and 
a DRAG component, yielding pulse shapes named Wah¬ 
Wah (Weak AnHarmonicity With Average Hamiltonian). 
Their results reveal that the WahWah ansatz works well 
if the targeted gate time is well-chosen. 

A key requirement for quantum computing on a 
large scale is for errors arising from leakage, address¬ 
ability, and frequency-crowding to fall well below a 
fault-tolerant threshold where quantum error correction 
(QEC) schemes [31, 32] improve with system size - as 
has been recently demonstrated by Kelly et al. [33] for a 
linear array of up to nine qubits. Gonservative estimates 
[34, 35] for this threshold state that the average error 
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per gate should be ;< 10“^, with lower errors also requir¬ 
ing far fewer physical qubits per logical qubit, and with 
leakage out of the computational subspace being an ad¬ 
ditional but not insurmountable impediment [36]. Thus, 
it is all the more pressing that these forms of error be 
suppressed. 

These problems become compounded when multiple 
gates need to be operated simultaneously. This operation 
is crucial to reducing computation times, which in turn 
allows for algorithm error linearly smaller in the number 
of qubits, as well as lower error budgets in QEC, due to 
error also accumulating during memory operations. In 
the presence of near-resonances in crowded systems on 
the order of 1% of the natural frequencies, pulses derived 
from adiabatic techniques are too short with respect to 
the decoherence time, whereas sideband modulation in 
the spirit of Ref. [28] disrupts concurrent gates with 
which it does not commute, leading to disruption of the 
desired final-time spectrum and even worse leakage than 
not including it. 

In this work, we rather consider decoupling the effect of 
the leakage transition using interference between the dif¬ 
ferent portions of the gate, being reminiscent of compos¬ 
ite pulses. This is accomplished by expanding the control 
pulses in terms of (Hanning) windows, which include a 
third harmonic that can be used to drive three sequen¬ 
tial rotations with alternating direction. The weights of 
the components are optimized numerically, yielding high- 
fidelity simultaneous X and/or Y gates with very short 
quantum speed limit, on the order of 27r/(5, where 6 is 
the detuning between the two most closely crowded fre¬ 
quencies. The waveform shapes found are similar to the 
CORPSE [37] pulse sequence, and share a robustness to 
deterministic and non-deterministic phase shifts during 
the pulse. 

This work is structured as follows: In Sec. II we intro¬ 
duce the system composed of two (3D) transmons being 
collectively driven through a cavity. We derive various 
frame transformations that allow for analytically quan¬ 
tifying the error. Within Sec. Ill we outline the analytic 
ansatz that is further optimized numerically and ana¬ 
lyzed with respect to certain properties, such as quantum 
speed limit and robustness, within Sections V and VII. 
In Sec. IV we outline the degrees of freedom and show 
samples of our pulses. An extended model for WahWah 
pulses is given in Sec. VI, achieving high fidelities for all 
gate times beyond a quantum speed limit. 


II. SYSTEM 


A. Model in the lab frame 

The computational subspace used for the qubit system 
in superconducting qubits is usually formed by the lowest 
two energy levels of an anharmonic oscillator with weak 
anharmonicity [16], although this choice is not manda¬ 
tory. In order to have a realistic model for the processes 


occurring in the system it is inevitable to consider at 
least the next higher energy level - which is referred to 
as a leakage level in the remainder of this work. Starting 
in the lab frame we can break down the Hamiltonian of 
our system, which is composed of two superconducting 
transmon qubits, into a (constant) drift part Hg and a 
controllable part He. Our bare two qubits are weakly 
anharmonic oscillators, described by Hq. Their energy 
level diagram is depicted schematically in Fig. 1. The lab 
Hamiltonian reads 


H = Ho + He 
2 
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This form of the Hamiltonian can be derived from quanti¬ 
zation of electrical circuits as shown in detail in Ref. [38] 
and in Ref. [39] for the two-level case. Given that the sys¬ 
tem is operated in the dispersive regime, i.e. cavity and 
qubit are sufficiently far detuned, the effective Hamilto¬ 
nian can be recast into the form of Eq. (1). Without loss 
of generality, we can assume our anharmonicities to be 
equal, i.e. Ai = A 2 = A. This will not lead to any de¬ 
generacies since the detunings of different levels remain 
non-degenerate and independently tunable. The perpet¬ 
ually used projection and generalized Pauli operators are 
defined as 


nf) = (2) 

= (3) 

3 

d;i'!^i = lj)(j-il('=) + b--i)(jf'=^ (4) 

= i\j){j - -i\j- (5) 

where the superscript (k) refers to either qubit 1 or 2. 
n{t) is a semiclassical dipolar interaction control with 
four quadratures total, paired into two sets, whereby each 
set’s carrier has frequency uJdi and uid 2 , respectively, so 
that 

^{t) = exi{t) cos {uJdit + (/i) -I- eyi{t) sin {ujdit + (/i) 

+ £x 2 it) cos {ujd 2 t + (l> 2 )+ £y 2 it) Sin {uJd 2 t + 4 > 2 ) ■ 

( 6 ) 

(k) 

The A) are associated to the relative strength at which 
the control field H(t) drives the leakage transition jl/j) o 
\2k) compared to the working transitions jOfc) <->■ jlfc). 
Numerical values for occurring frequencies and A^^^ are 

(k) . '' 

given in Tab. I whereby the A) are approximated to be 
•y/J since the anharmonicities are assumed to be small 
(hence eigenstates are close to those of an harmonic os¬ 
cillator) and since the transmons are far detuned from the 
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Qubit 1 Qubit 2 

FIG. 1. Energy level diagram of the two qubits. The leakage 
transition II 2 ) |22) of qubit 2 is only slightly detuned by 5, 
which we call the crowding frequency in the remainder of this 
work, from the working transition |0i) |li) of qubit 1, i.e. 
nji + (5 = 0^2 + A. The dashed lines indicate the level structure 
if the transmons where harmonic oscillators with frequencies 
uj\, UJ 2 - We collectively drive the system with pulses at ujdi 
and ujd 2 - 


driving resonator [24]. In principal, the relative phases 
4> i ,2 between envelope and carrier need to be taken into 
account. However, since we apply the rotating wave ap¬ 
proximation (RWA) throughout following sections, those 
phases become irrelevant [24, 40]. 


(k) 

Particular choices of can lead to frame representa¬ 
tions with different meaning and application purposes, 
such as 

= ju>d , rotating frame with ujd 

+ IsSp , interaction picture. 

Transforming to a different frame along the lines of Eq. (8) 
will also affect the unitary time evolution operator ac¬ 
cording to 

U^(t) =R(t)U(t)Rt(0), (9) 

where U^(t) is the time evolution operator inside the 
rotating frame. 


1. Rotating frame 

In order to have a better insight in the dependence 
of our system on the two driving frequencies we choose 
to look at the rotating frame with ' = juJdi- For 
more controllability we pick our carrier frequencies with 
arbitrary fixed detuning Ai 2 from the qubits’ resonance 
frequencies: 


Wdl — W1 -b Ai (10) 

Wd2 = W 2 -b A 2 (11) 


TABLE I. Definition of the system parameters’^. 



Qubit 1 

Qubit 2 

Unit 

oJkl'in 

5.508 

5.903 

GHz 

A/27r 

-350 

-350 

MHz 

& 1‘1'K 


45 

MHz 


1 

1 


\(k) 

^2 

V2 

V2 



^ Taken from Ref. [28]. Suggested for experiments by the group of 
Leo DiCarlo at TU Delft. 


B. Frame transformations 


As Hq and R are diagonal, our new Hamiltonian reads 

H« = Ho -b R^HcR + (12) 

In the remainder of this work we will use the shorthand 
notation 

= exj{t) + ieyj{t) and (13) 

x(t)=Hi(t) + e*^‘H2(t), (14) 

with 7 = ujdi — U!d 2 = A — d -b Ai — A 2 being the detun¬ 
ing between both carrier tones. Performing the RWA, 
i.e. neglecting terms that oscillate with zL 2 uJdi ,2 and 
± {uJdi + a^d 2 )j leads to the rotating frame Hamiltonian 


Since we are interested in controllability of our system 
by microwave fields (amplitude and carrier) we choose 
to work in different frames, such as a frame rotating at a 
specific frequency in order to remove any intrinsic oscilla¬ 
tion from the controls. Such a transformation to a frame 
|$fl) = R'f 1$) rotating at frequencies is achieved 

by applying 


H-^ = -Aift^^^ -b (A - 2Ai)n^^^ 

-b (d - A - Ai)nf ^ -b (2(5 - A - 2Ai)n^^^ 


m = 




( 7 ) 


-bh.c. 


(15) 


according to the transformation rule 

H« = Rl'HR-biR^R. (8) 


Evidence for the validity of the RWA is that the system 
frequencies in the rotating frame are of the order of A, 
thus more than an order of magnitude smaller than the 
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qubit frequencies uji^ 2 - At this point we want to briefly 
comment on why a simple DRAG solution alone will 
in general not improve results. The main idea behind 
DRAG is to move to an adiabatic frame via a transfor¬ 
mation V, satisfying V(0) = V(tg) = 1. For simplicity 
we will restrict to only one driving field, hence equiva¬ 
lent to fi 2 (t) = 0 Jd 2 = A 2 = 0 in Eq.(15). Essentially, the 
transformation matrix will then be 


V (t) = exp 




A 

n 

j=i 


A A:- 'cr 


(2)^y(2) 

0 1 


(16) 


computational subspace) is given by 


Ured(tg) = 6*“^ exp 




(18) 


with an arbitrary, yet unimportant global phase 4>. A 
complex rotation angle 9 gives rise to arbitrary rota¬ 
tions around the X/Y axis. The real part Re{0} will 
rotate around the X-axis whereas the imaginary part 
lm{0} rotates around the Y-axis. We will develop and 
design pulses starting from the interaction frame, so that 
throughout Sec. IIIB the final time evolution operator of 
the two transmon system is approximately given by 


A key to sufficiently simple DRAG solutions is fast con¬ 
vergence of a series expansion with respect to an expan¬ 
sion parameter rj = ExiIk 1 [23]. Performing the cor¬ 
responding adiabatic expansion to first order will lead 
to multiple possible values for k and a compensation 
quadrature e^i oc {1/K)exi- Eor instance, k = d removes 
the leakage from qubit 2. However, this choice will leave 
in considerable phase error on both qubits and slightly 
increase errors from the other (less) crowded frequencies. 
In principle, similar problems occur for other choices of k 
so that a simple DRAG ansatz will not be able to effec¬ 
tively address frequency-crowding issues. One can also 
optimize within a set of higher derivative functions, as 
was done in Ref. [25]. However, we will show that a 
simpler non-adiabatic basis will exist for our problem. 



1 

2 


/O 


o\l 


0 

0 

VO 

0 

0/J 


(19) 


rotating qubit j about an angle 9j and neglecting phases 
on the leakage level. Gonsequently the gate we actually 
implement in the lab frame is obtained by utilizing Eq. (9) 
to transform Eq.(19) back to the lab frame. Eventually 
this will lead to Z errors (relative phase shifts) that are 
quantified in Sec. VHI together with techniques to com¬ 
pensate them. 


B. Magnus expansion 


2. Interaction frame 


1. Basic idea 


Transforming the lab frame Hamiltonian according to 
Eq.(8) to the interaction frame yields the interaction 
Hamiltonian 


An arbitrary Hamiltonian H(t) induces a time evolu¬ 
tion after a time tg according to the time evolution op¬ 
erator 


1 = 1 

-I- h.c. 


(17) 


where we have introduced the eigenfrequencies = 

-Ai, 4^^ = A-Ai,(5f^ =(5-A-Ai and 4^^ =J-Ai. 
In the remainder of this paper we will be working in the 
frame defined by Eq. (17). 


III. SIMULTANEOUS SINGLE-QUBIT GATES 


U(tg) = T exp 



( 20 ) 


where T is the time-ordering operator. The latter ac¬ 
counts for the general impossibility to compute a closed 
analytical form for U(t) due to the fact that H(t) in 
general does not commute with itself at different times. 
An approximation for the final unitary U(tg) can be ex¬ 
pressed in terms of the Magnus expansion [41, 42] 

U(tg) = exp Okitg)^ ■ (21) 


A. Target evolutions 

At first, we focus on implementing single qubit rota¬ 
tions about an rotation angle 0 € C so that in the lab 
frame the desired unitary for each qubit (reduced to the 


The ©fe’s are hermitian matrices generated by nested 
time integrals over nested commutators of the Hamilto¬ 
nian at different times. A major advantage of the Mag¬ 
nus expansion is the straightforward treatment of time¬ 
ordering. The first two orders of the expansion are given 
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by 


©0(^3) 


©1(^3) 



( 22 ) 


Convergence of the series is not always guaranteed. For a 
differential equation Y{t) = A{t)Y{t) in a Hilbert space 
7i with boundary condition F(0) = 1 it is proven [42] 
for bounded operators A{t) that /p dt ||^(t)|| < tt is a 
sufficient condition to guarantee convergence in the time 
interval [0,T). Nevertheless, this condition is only a suf¬ 
ficient one so that the series may still converge for t > T 
even if the previous criterion is not satisfied. In case of 
Y(t) being a normal operator, in particular a unitary one, 
the series also converges for infinite dimensional prob¬ 
lems. 


respectively. In the above equations the correspond¬ 
ing complex conjugate versions will also hold. For fol¬ 
lowing calculations we introduce complex-valued ampli¬ 
tudes ai _2 so that x(0 = aifii(t) -I -02 with 

rescaled pulses Hi^ 2 - These rescaled pulses are further 
decomposed into their real and imaginary parts, i.e. 

= ixjit) + From Eq.(17) it follows that 

Eqs.(24a) and (24c) are the working transitions of both 
qubits, whereas the other two equations belong to their 
leakage transitions. Semiclassically speaking, conditions 
(24b) and (24d) state that there must be no spectral 
power at the unwanted transitions whereas conditions 
(24a) and (24c) are instances of the area theorem. Terms 
arising from © 1 (^ 3 ) will mostly contribute to diagonal el¬ 
ements, thus characterizing Z errors in the final unitary. 
Those Z errors mainly stem from level shifts and popu¬ 
lation leaking out of the computational subspace while 
driving the system. 


C. Hanning windows 


2. Lowest order conditions 

Requiring the lowest order Qo{tg) of the Magnus ex¬ 
pansion to implement the gate of Eq.(18), i.e. 

tjpitg) = exp (^-i0o(tg)j (23) 

gives a first idea about pulse characteristics. That way, 
errors may arise from higher orders in the expansion 
which can partially be related to an AC stark shift or 
a Bloch-Siegert shift [43, 44]. Such errors can for in¬ 
stance be reduced by either requiring higher orders to 
vanish or to lead to global phase factors. Since the in¬ 
teraction Hamiltonian can be written as a direct sum 
H/ = Hj^^ 0 Hj^^ and using the identity exp(A 0 H) = 
exp(A) 0 exp(i3) for A = H^^^ and B = H©^ (superscript 
(k) refers to qubit k) we find the conditions 

^9 

J dt = di (24a) 

0 
^9 

J dt = 0 (24b) 

0 
^9 

Af ^ J dt = 02 (24c) 

0 

^9 

A^^^ J dt = 0 (24d) 

0 

to implement unitaries having the form of Eq.(I9). Here, 
01,2 G C are the rotation angles of qubit I and two. 


Gaussian pulse shapes are widely used to implement 
single-qubit rotations since they are smooth and have a 
limited bandwidth in excitation spectrum. Taking into 
account leakage levels and spectral crowding issues, new 
analytic methods such as DRAG and WahWah have been 
developed to counteract those challenges. However, in 
case of simultaneous gates, the large number of con¬ 
straints makes the choice of basis functions (e.g. Her- 
mite, Gaussian, Error function) especially important for 
considerations such as smoothness, boundary conditions 
and number of parameters. 

Yet, WahWah pulses only suppress leakage at the end 
of a gate, so that simply applying a WahWah control 
resonant with each transmon’s logical transition will not 
allow simultaneous gates since there will be too much 
crosstalk between both physical qubits. We will choose 
to use the Hanning window functions [45] 


N 

n—1 


1 — COS 



(25) 


which guarantees smoothness and also satisfies the 
boundary conditions for the pulse (start/end at zero). 
Utilizing this family of functions, we will also show that 
small errors at short times using only a small number of 
parameters can be obtained. 

Hanning windows are expected to perform well in large 
part due to having both features that include smooth¬ 
ness (which helps to enforce adiabaticity) and composite- 
pulse structure that allows cancellation of errors between 
different (e.g. diabatic) parts of the pulse. Moreover, 
the composite-pulse nature (very similar to the compos¬ 
ite decoupling sequence, CORPSE [37]) can have intrin¬ 
sic robustness to frequency miscalibrations, noise and 
time-dependent shifts (e.g. the Stark shift caused by 
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FIG. 2. Control amplitudes for 26ns (a) and 42ns (b) together with their respective spectra (c,d). For longer gate times some 
controls resemble simple Gaussians, whereas they show characteristics of composite pulses especially at shorter gate times to 
account for large Z errors during the sequence. The spectra indicate that Eqs.(24) are essential to implement the desired 
rotations in crowded systems. 


off-resonant drive of the frequency crowded level), see 
Sec. VII. 

The performance of these pulses will be discussed and 
compared to Gaussians, derivative-based functions and 
the WahWah ansatz in the next sections. 


1. Constraints 


From the lowest order Magnus expansion we find that 
the finite Fourier transform 


s[n,p) 


dt Ct{t) 


jpt 


(26) 


gives - to lowest order - a good intuition about the pulses’ 
characteristics. Solving Eqs. (24a) and (24c) for the 
amplitudes so that the controls implement 0i-rotations 


yields 


Oi = 


02 = 


'-025(^2, a-^-A a) 


' jh _ 01 A- ki) 

5(122,-Ai) 


X 5(122,-A 2 )- 


5(f2i, (5 - A - Ai)5(f22, A - <5 - A 2 ) 


-1 


5(f2i,-Ai) 


(27) 


Since the two logical transitions are detuned by 5 — A ^ 
400 MHz from each other, it is valid to neglect the in¬ 
fluence of both controls to the working transition they 
are not near-resonant with if their spectrum shows suf¬ 
ficient decay away from resonances, i.e. if for instance 
5(i2i,(5 — A — Ai)/5(f22, — A 2 ) 1. Thus, we simply 

want each control to have an area complying with the 
desired angle of rotation 0^, hence 

dt f2i(t)e-*^’* j . (28) 

In fact, we observe that the approximation giving rise 
to Eq. (28) leads to the same performance as the exact 
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amplitude conditions in Eqs. (27). Finally, after picking 
an ansatz for the tilded x/y-control shapes, the controls 
in Eq.(6) are obtained as 

Sxiit) = Re |ai , eyi(t) =Im|aifii(t)| 

£x2it) = Re|a2^2(t)} , £y2{t) = Im|a2fi2(0} ■ 

(29) 


IV. PULSE SHAPING 


A. Degrees of freedom & Optimization 


To quantify our results we make use of the common 
overlap gate fidelity 



Tr 




t 

target 



2 


(30) 


which is insensitive to global phases. Here, the trace 
is taken over the compound computational subspace 
{|00), |01), |10), |11)} with d = 4 being its dimension. 

Investigating the performance of multiple ansatzes for 
the tilded controls reveals that it is satisfactory to choose 
rixj to be a superposition of the first three Hanning win¬ 
dows, i.e. setting iV = 3 in Eq. (25). Since one out of the 
three coefficients per control can be incorporated into the 
corresponding amplitude Ui, there are effectively four co¬ 
efficients plus two detunings Ai _2 that determine the full 
set of parameters. To account for higher order errors, 
such as level shifts and Z errors, we use the Nelder-Mead 
(NM) simplex algorithm [46] to obtain (locally) optimal 
values for the degrees of freedom. Since the tendency 
to get trapped in local extrema increases with the num¬ 
ber of parameters [47] it is at least necessary to sample 
over various initial conditions and pick the best found 
solution. 

We can expect that generalization to more than two 
qubits if additional qubits do not add crowding on the 
order of 6 is straightforward. In this case, all waveforms 
can be supplemented by DRAG to reduce the impact of 
crosstalk to other qubits - scaling linearly with the num¬ 
ber of qubits [25]. Even in case of multiple frequency- 
crowded qubits, one can in principle utilize higher order 
Hanning windows and/or additional quadratures to still 
achieve precise simultaneous gates. Nevertheless, this 
will ultimately increase the parameter space and thereby 
hamper the optimization. There are indications [48] that 
one might still observe reasonable convergence if operat¬ 
ing far enough away from the quantum speed limit. A 
detailed study of this aspect is an area of future research. 


B. Sample pulses 

Examples of typical control shapes and their finite 
Fourier transforms are shown in Fig. 2. For shorter gate 


-kZ 



-z 


FIG. 3. Visualization of the state evolution \tp 2 } with 
|V’2(0)) = IO 2 ) for the 26ns pulses shown in Fig. 2(a). Due 
to the composite-like structure Z errors are corrected by ap¬ 
proximately applying a -|-27r/3 —> —tt/S —>■ -|-27r/3 pulse. The 
color encodes the amount of population inside the computa¬ 
tional subspace. Red means all population is inside {|0), |1)} 
space, while blue indicates leakage by about 25%. 


times, i.e. 26ns in part (a), the Sx 2 ,y 2 controls have 
magnitudes much larger than those required for a sim¬ 
ple TT rotation. That is, these amplitudes push the pulses 
into the diabatic regime, with large intermediary leakage. 
Moreover, they show a composite-like structure very sim¬ 
ilar to a well-known method for robust pulses (CORPSE) 
[37], approximately steering a -|-27r/3 —)■ —7r/3 —)■ -|-27r/3 
evolution on qubit 2, which is illustrated on the Bloch 
sphere in Fig. 3. We will show in Sec. VH that this com¬ 
posite structure enables robustness to frequency miscal- 
ibrations, noise and time-dependent shifts. 

Our analysis shows that this composite-like shape, 
which is especially important for shorter gate times, is 
needed to achieve gate errors sufficiently small to make 
the pulses feasible for experimental implementations. 
Since leakage out of the computational subspace is not 
suppressed during the drive, there must be a precise inter¬ 
play between exi,yi and 6 x 2 ,y 2 that ensures the correct ra¬ 
tio between how much population is pulled back into/out 
of the computational subspace by 6 x 2 ,y 2 while exi,yi drive 
a rotation on qubit 1 (still affecting the leakage transition 
on qubit 2 due to spectral crowding). Owing to severe 
leakage out of the {IO 2 ), II 2 )} subspace and due to level 
shifts, there will be relative phase errors, manifesting as 
rotations around Z axis. This kind of error can be com¬ 
pensated by properly adjusting the interplay between all 
controls with the help of the NM algorithm. For longer 
gate times, such as 42ns in Fig. 2(b), the composite struc¬ 
ture is not that essential anymore: in this domain the 
type of control solutions becomes different. While exi,yi 





B. Simultaneous rotations 


Gate time tg (units of 2n/S) 



FIG. 4. Final leakage error of qubit 2 at different gate times 
for Hanning controls. Except for some gate times we achieve 
leakage errors below 10“® after ~ 25ns. The non-monotone 
shape of the curve is most likely due to the system being 
overconstrained. 


resemble simple Gaussian controls with maximal magni¬ 
tude at t = tg/2, their counterparts £x 2 ,y 2 are dominant 
at the beginning and end of the pulse. This also illus¬ 
trates the interplay between all controls, leading to the 
correct cycling in and out of leakage subspace of qubit 
2. The spectra in Fig. 2 illustrate that the lowest order 
conditions in Eqs.(24) are essential to steer simultane¬ 
ous rotations. Nevertheless, spectral arguments are not 
sufficient because of non-negligible higher order errors. 


In Fig. 5 we show performance of resonant (dashed ma¬ 
genta) and off-resonant (dash-dotted blue) Hanning con¬ 
trols together with a naive Gaussian approach, one per 
qubit (solid black), and a derivative-based ansatz (dot¬ 
ted red). Due to the finite bandwidth of the controls, 
at least the control resonant to uji will have severe spec¬ 
tral power at the leakage transition 6 of the other qutrit, 
thus populating its |2) state and consequently leading 
to high gate errors in general. For instance, the gate 
error of Gaussian pulses slowly decreases exponentially 
with gate time. Derivative-based pulses perform roughly 
one order of magnitude worse compared to resonant Han¬ 
ning controls, which is due to less control over phases via 
composite-like shapes. Clearly, using at least the first 
three Hanning windows {N = 3) is sufficient for precise 
pulses in a frequency crowded system. Resonant controls 
enable average gate errors between 10“^ — 10“^ while de¬ 
tuning the carrier slightly (by a few kHz to MHz) from 
the qubits’ resonances further decreases the error by two 
orders of magnitude since there is considerably more con¬ 
trol over level shifting effects (Stark, Bloch-Siegert) than 
in the resonant case. In fact, it may also be sufficient 
to apply only the controls Sx 2 ,y 2 off-resonantly and leave 
Sxi,yi resonant with wi: this will slightly increase the av¬ 
erage gate error, however still being roughly 1.5 orders of 
magnitude better than having all controls on resonance. 
Nevertheless, the gain in fidelity by introducing more de¬ 
grees of freedom comes at the cost of slower convergence, 
as has been stated in section Sec. IV A. Thus, fewer con¬ 
trols may be preferable in practical situations. Fig. 6 il¬ 
lustrates that our ansatz allows precise implementation 
of all simultaneous rotations around the X/Y axes. We 


V. QUANTUM SPEED LIMIT 
A. Leakage error 

At first, we want to demonstrate that our method ef¬ 
ficiently suppresses leakage arising from spectral crowd¬ 
ing in our system. Fig. 4 illustrates that we are able to 
achieve leakage errors (due to the critical |22) state) that 
are well below 10“^ after a gate time ~ 25ns, indicat¬ 
ing the relationship to the speed limit found in Sec. V C. 
The observed limit of approximately 25ns for sufficiently 
suppressing leakage is largely explained by a bandwidth 
argument: For pulses of duration tg < 27r/(5 ~ 25ns it 
is expected from spectroscopy that transitions with fre¬ 
quency close to 6 will be excited, hence notable leakage 
out of the computational subspace occurs. Besides its 
impact on gate errors, leakage plays an important role in 
quantum error correction schemes. The majority of those 
assume that there will be no information leakage out of 
the computational subspace [36], making them vulnera¬ 
ble to leakage errors. 


Gate time tg (units of 2nj5) 


0.5 0.75 1 1.25 1.5 1.75 2 2.25 



FIG. 5. Performance of Gaussian pulses (each resonant to 
^ 1 , 2 ) , optimized Hanning-based pulses and derivative-based 
controls. Off-resonant controls yield an improvement by 2-3 
orders of magnitude compared to their resonant counterparts. 
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Gate time tg (units of 2tt/6) 



Gate time tg (ns) 


FIG. 6. Average gate error as a function of gate time tg for 
various simultaneous rotations that are important for appli¬ 
cations, for instance, in the AllXY sequence [49]. 


Logio(l - #) 



Detuning (5/27r (MHz) 

FIG. 7. Average gate error as a function of gate time tg 
and crowding frequency S. The solid black line indicates the 
best fit for the speed limit (referring to 10“"^ error) which is 
essentially t™‘“((5) oc 1/5. This scaling indeed complies with 
the limit up to which leakage can be sufficiently be suppressed, 
as explained in Sec. V A. 


focus on rotations that are part of the AllXY sequence 
[49] which can be used to tune up quantum systems and 
identify various error sources in order to underline the 
close connection between our results and recent experi¬ 
mental requests. Simultaneous gates involving identities 
are also easily implemented using Hanning-based pulses. 
However, in section Sec. VI we present a generalization 
of Ref. [28] that is exactly designed for this purpose. 
Various other single qubit gates, such as the Hadamard 
gate 


VI. SEPARATE SINGLE-QUBIT ROTATIONS 

As mentioned in the previous section, our Hanning- 
based pulses can also be utilized to implement 1 gates. 
Nevertheless, the WahWah technique [28] was designed 
for that purpose. There, only two control quadratures 

e^{t)=A„e ^ |l - cos I (32a) 

£y(.t) = (32b) 


tJn 





(31) 


can be decomposed into products of single qubit rota¬ 
tions. For instance, Uh = V^Vr /2 (up to a global phase) 
may be implemented in roughly 40ns - even in case of 
very strong crowding which is primarily analyzed in this 
work. 


C. Dependence on spectral detuning 


Because of the inevitable crosstalk between both trans- 
mon qubits, one expects a gate time (as a function of 
crowding S) at which precise control (we refer to the 10“^ 
error threshold) is no longer possible. Since excitations of 
transitions with frequency S approximately scale oc 1/5, 
we expect the quantum speed limit to scale in the same 
manner. In Fig. 7 we show the optimal average gate error 
as a function of gate time tg and crowding frequency S. 
The black line depicts the best fit, substantiating that 
the speed limit t™“(5) ~ 2n/6. 


resonant with wi are exploited to drive a 7r-rotation on 
qubit 1 while leaving alone the second one, i.e. U(tg) = 
X (g) 1. To this end, a Gaussian is sideband-modulated 
with frequency Wx that allows cancellation of leakage er¬ 
rors from qubit 2. The standard deviation is chosen as 
a = tg/6 to let the pulses smoothly start and end close 
to zero. In order to prevent qubit 1 from leaking, the in- 
phase control Ex is supplemented with DRAG whereby 
the factor of 2 in the denominator stems from the ab¬ 
sence of phase control [24]. 

As pointed out in Ref. [28], there will be a relative 
phase error on qubit 2, eventually increasing the gate 
error measured via the fidelity in Eq.(30). Therefore, the 
fidelity functions 


= 7 


®avg - 2 


Tr{|0,i).|l,i>} {uJ;.U(tg)| 

(‘J’p.O) + ‘J’l*,!)) 


2 


(33) 

(34) 


are introduced. Here, the partial trace is taken over 
all states where the second qubit is either in |0) or |1). 
Hence, a high outcome of 4>|, implies that qubit 2 starts 
and ends in state |i). Owing to this, the averaged reduced 
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Gate time tg (units of ‘I'KjS) 

0.5 0.75 1 1.25 1.5 1.75 



Gate time tg (ns) 


FIG. 8. Gate error for numerically optimized WahWah pulses 
and pulses defined by the model in Eq.(35), compared to a 
Gaussian and the original WahWah pulses. We observe a 
speed limit for the optimized pulses at ~ 0.75 2tx ji. 


means to account for such imprecise knowledge is the Ad- 
HOC protocol [50] which combines open- and closed-loop 
optimal control, enabling enhancements in gate fidelities 
by an order of magnitude. However, since errors can in 
general occur during/in between operations of the hard¬ 
ware, e.g. 1// noise, it is desired to work with pulses 
that show intrinsic robustness against certain errors and 
thereby go beyond closed-loop control methods. 


A. Filtering effects 

When applying optimal control pulses in an experi¬ 
ment, their shape will be altered according to the transfer 
function of used hardware. For instance filtering effects 
due to finite bandwidth of waveform generators is a cru¬ 
cial part of those transfer functions. By transforming the 
controls in Eq.(29) according to 

+ 00 +00 

Mfiit.(i) = ^ J dt' J duj E(w)e*(‘“*')‘^'Uunfiit.(t') (36) 


fidelity $avg is a way to estimate performance insensitive 
to relative phases in qubit 2. Originally, a modulation 
with Wj; = S/2 was suggested and turned out to be suc¬ 
cessful if the gate time was well-chosen. From numer¬ 
ical studies, we derive the optimal sideband frequency 
= ^x/S in terms of tg = tgS/2Tr as 

^x{ig) = 

2.3 erf h.ld^ig - |) , f < ig < f 

2-3 erf (^)+0.41 (4-1), ig>l 

(35) 

The model works well in a large range of (5/A and incor¬ 
porates a speed limit for the ansatz in Eqs.(32), which is 
approximately given by 0.75 2Tr/6. For gate times close to 
the speed limit, there is a nonlinear dependence propor¬ 
tional to the error function erf (x). After passing through 
the nonlinear region tg > 5/4 the sideband frequency 
scales linearly with the gate time. In Fig. 8 we illustrate 
performance of fully optimized pulses (solid black line), 
the piecewise model in Eq.(35) (dash-dotted blue line), 
the pure linear model for (dashed magenta line) and 
a naive Gaussian approach (dotted gray line). 

Phase errors in qubit 2 can be tracked and corrected 
by methods presented in Sec. VIII. 

VII. ROBUSTNESS ANALYSIS 

An important property for experimental application of 
control pulses is their robustness against mismatches in 
for instance carrier frequencies or pulse envelopes, as well 
as filtering effects due to experimental hardware. One 


we model hardware filters. Here, F(u}) is the response 
function of the filter which is assumed to be Gaussian, 
i.e. F{u}) = exp(— w^/clIq) whereby ojq/2tt = 425.4 MHz 
(approximation for Tektronix AWG50I4 [40]). Eq.(36) 
describes a non-causal filter, which was found [40] to be 
a better approximation than a causal one. Applying this 
filter without further corrections has - up to some ex¬ 
ceptions - no tremendous effect on fidelity, as shown in 
Fig. 9. We find that fine-tuning the control amplitudes’ 
magnitudes almost corrects for all additional error, either 
by precomputing the effect of the filter, or making use of 
closed-loop (Ad-HOC) type optimization in experiment. 
In general, if the transfer function of experimental hard¬ 
ware is known, its inverse can be incorporated into the 
input controls to achieve better performances [40]. 

B. System parameters 

Besides filtering, also uncertainties in system parame¬ 
ters such as qubit frequencies may give rise to additional 
errors. The two most important frequencies being cru¬ 
cial for the system dynamics are the anharmonicity A 
and most notably the crowding frequency 5. It is ex¬ 
pected that the optimal solutions will strongly depend on 
those two frequencies, giving rise to the question about 
robustness of the solutions against imperfect knowledge 
of 5 and A. Fig. 10 depicts the average gate error un¬ 
der variations of A and 5, showing that the solutions 
are very robust against mismatch in anharmonicity and 
pretty robust against uncertainties in 5: Deviations of <5 
by ±1.5% and ±4% can still lead to gate errors of 10“^ 
and 10“^, respectively. The observed asymmetry in A 
stems from the fact that A < 0. A negative deviation of 
A means that |A| decreases, so that leakage errors inside 








11 


Gate time tg (units of 2-11/5) 



FIG. 9. Performance of optimal control shapes under Gaus¬ 
sian filtering compared to the unfiltered case. Except from 
some gate times, almost all increase of error can be compen¬ 
sated by adjusting the pulse’s magnitudes. 


Logiofl - #) 



Deviation A/ 27 r (%) 


FIG. 10. Error landscape for a X (g) X gate at 30ns. The 
landscape illustrates that our pulses are very robust against 
deviations of the anharmonicity A. Although 5 is expected to 
be the most crucial parameter, we are able to achieve errors 
< 10“^ if 5 is known up to only ±4%. 


qubit 1 become more likely whereas a positive deviation 
increases the distance between |1) o |2). Note that in 
general, also timing jitter from the arbitrary waveform 
generator will affect system performence. As it will be 
small and can by rescaling usually be compared to sys¬ 
tem parameter uncertainty, we expect its impact to be 
similarly small. 


VIII. CORRECTING RELATIVE PHASES 


As mentioned before (see Sec. Ill A) we introduce Z 
errors and global phases in our lab frame gates by the 
ansatz we choose. Since global phases are irrelevant we 


do not have to account for them; however, Z errors should 
be corrected. We will only focus on the correction inside 
the computational subspace because phases on the leak¬ 
age levels are unimportant. Neglecting global phases, we 
calculate the Z errors we make from Eq.(9) as 


exp 


. (Wdl — Al)tg 
* 2 


® exp 


. (wdi -h (5 — A 
* 



This can easily be corrected if a Z control is available on 
both qubits by applying controls Zx and Z^ for a duration 
T so that their areas are 


l\t Zx{t) = 

Tdt z,{t) = 

./n ^ 


(38) 

(39) 


Alternatively it is also possible to account for this error 
by properly adjusting the phase of the preceding gate. 
This procedure is similar to what is named phase ramping 
[24, 51] or frame compensation [25]. Essentially the XY- 
plane is rotated by an appropriate angle to account for 
the error made. 


IX. CONCLUSIONS 

We have investigated independent control of two trans- 
mon qubits coupled to the same cavity. Typically, high 
fidelities in a system like this are hindered by spectral 
crowding, where for instance a harmful leakage transi¬ 
tion is close in frequency to a logical transition (detuning 
S) and thereby renders individual addressability impos¬ 
sible. Our ansatz, based on a superposition of Hanning 
windows, achieves gate errors well below 10“^ while be¬ 
ing subject to a quantum speed limit slightly below 2Tr/6. 
We have shown that detuning the controls from resonance 
by a few kHz to MHz can again decrease the gate error 
by two orders of magnitude. 

Our pulses outperform adiabatic methods and Wah- 
Wah pulses for fast and precise implementations of si¬ 
multaneous single-qubit rotations. Moreover, the ansatz 
theoretically allows for an easy generalization to more 
than 2 qubits, especially if additional ones do not add 
crowding on the order of 6. We have adressed possible 
limitations in Sec. IV A. 

Additionally, we have presented a model that general¬ 
izes WahWah pulses to arbitrary gate times and crowd¬ 
ing frequencies, as long as not aiming at times below the 
quantum speed limit ~ 0.75 2Tr/6. 

We have shown that pulse oscillations occur on 
timescales that are unsusceptible to filtering effects and 
that their composite-like shape gives rise to robustness 
against off-resonance errors, such as mismatches in char¬ 
acteristic frequencies. Moreover, the pulses are analytical 
and allow for easier debugging and benchmarking com¬ 
pared to fully numeric control shapes. 
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